Critical phenomena of the hard-sphere lattice 
gas on the simple cubic lattice 

Atsushi Yamagata 

Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro- 

ku, Tokyo 152, Japan 

and 

Department of Physics, Tokyo Metropolitan University, Minami-ohsawa 1-1, 
Hachioji-shi, Tokyo 192-03, Japan 

Running title Critical phenomena of hard-sphere lattice gas 

Keywords Hard-sphere lattice gas, Critical phenomena, Monte Carlo method 

PACS classification codes 02.70.Lq, 64.60.Cn, 68.35.Rh 

Abstract 

We study the critical phenomena of the hard-sphere lattice gas on the simple 
cubic lattice with nearest neighbour exclusion by the Monte Carlo method. 
We get the critical exponents, = 0.313(9) and 7/z/ = 2.37(2), where /? 
is the critical exponent for the staggered density, 7 for the staggered com- 
pressibility, and V for the correlation length. 



1 Introduction 



Atoms of a hard-sphere lattice gas occupy sites of a lattice and interact with 
infinite repulsion of nearest neighbour pairs. The grand partition function is 
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N 



MN), 
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where z is an activity and Zv{N) is the number of configurations in which 
there are N atoms in the lattice of V sites. At 2; = +00 a ground state 
configuration is that the atoms occupy all the sites of one sublattice and the 
other is vacant. There is no atom at z = 0. A continuous phase transition 
occurs at a critical activity. There are many studies of the system by various 
methods: series expansions ||, |], finite-size scaling and transfer 

matrix 0, Bethe and ring approximations [0, transfer matrix ||^, ^, |T^], 
corner transfer matrix and series expansions [jll|, exact calculations |T^, and 
Monte Carlo simulations fl^, 

The critical activity is obtained on various lattice. The square lattice: 
3.80(2) 0], 3.80(4) g, 3.7966(3) fig, 3.7962(1) and 3.796255174(3) 
The triangular lattice: 11.12(11) [|], 11.05(15) [|], and ^(11 + 5v4) = 
11.090... The simple cubic lattice: 1.09(7) § and 1.0588(3) H. The 
body-centred cubic lattice: 0.77(5) P and 0.73(2) |T6|, [1^. The critical ex- 
ponents on the square lattice are as follows. The thermal exponent, yt, is 
1.000000(1), the magnetic exponent, yn, is 1.875000(1) H, f3/u = 0.125(5), 
and '-f/u = 1.74(2) ||13| where (3, 7, and u are the critical exponents for the 



order parameter, the staggered compressibility, and the correlation length, re- 
spectively. There are exact values on the triangular lattice [ll2| , Q: a = 1/3 
for the compressibility, f3 = 1/9, fi = 5/6 for the interfacial tension, and 
u = 5/6. The value, 0.61071(2) pi, of the fourth-order cumulant pi of the 



order parameter on the square lattice agrees with those of the Ising model, 
0.611(1) 0, 0.6103(7) 



|2g], 0.611(1), and 0.610690(1) ]I9[. These results 
support that the hard-sphere lattice gas on the square lattice belongs to the 
universality class of the two-dimensional Ising model and that on the triangu- 



lar lattice does to one of the three-state Potts model in two dimensions 
There is no result of critical exponents in three-dimensional hard-sphere lat- 
tice gas. 

We carry out Monte Carlo simulations of the hard-sphere lattice gas on 
the simple cubic lattice. Using the finite-size scaling, we estimate the critical 



1 



exponents. In the next section we define physical quantities measured. In 
section || we describe procedure of analyses and present Monte Carlo results. 
A summary is given in section ^. 



2 Monte Carlo simulations 

We use the Metropolis Monte Carlo technique ^ to simulate the hard- 



sphere lattice gas (P on the simple cubic lattice of V sites, where V = L x 
LxL (L = 2xn, n = 2,3, ... , 30), under fully periodic boundary conditions. 
According to Meirovitch |T3[, we adopt the grand canonical ensemble. The 



algorithm is described in the references [13, 14 



We start each simulation with a ground state configuration at the critical 



activity, Zc = 1.0588 [|^. The pseudorandom numbers are generated by the 
Tausworthe method We measure physical quantities over 1.2 x 10^ 

Monte Carlo steps per site (MCS/site) for L > 44 and 10^ MCS/site for 
L < 42 after discarding 5 x 10'^ MCS/site to attain equilibrium. We have 
checked that simulations from the ground state configuration and no atom 
one gave consistent results. Each run is divided into M blocks. Let us the 
average of a physical quantity, O, in each block {0)i; i = 1,2, ... , M. The 
expectation value is 

1 



The standard deviation is 



1=1 



1/2 



A{0) = ({O)^ - {O) /VaT^. 



We have taken M = 12 for L > 44 and M = 10 for L < 42. 
Let us define an order parameter by 

R = 2{Na-Nb)/V 

where A^a (^b) is the number of the atoms in the A (B)-sublattice and 
A^ = Aa + Ab- We measure the staggered densities. 



m 
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and 



mt' = (i?2)i/2, (3) 
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the staggered compressibilities, 



= V{{R^) - {\R\)^)/4 (4) 

and 

X^' = V^(^/4, (5) 
and the fourth-order cumulant of R |]20| , 

U = 1-^W)/{W- (6) 

3 Monte Carlo results 

We estimate a critical exponent and an amplitude by using the finite-size 



scaling [pSl , P9| , |30|| . For a physical quantitiy, O, we use the nonlinear chi- 



square fitting fSlI] with a function of L, 

0{L) = ALP, 

where p and A are fitting parameters, p = —(3/v when O = or m)' . 
p = ^jv when O = oi' X^' ■ We calculate the chi-square per degrees 
of freedom, x^/DOF, and the goodness of fit, Q, 0| for the data set of the 



sizes h — ivjnin; -^min + 2, . . . , Ljiiax~2, and -Z^max- The values of Ljiiin S-Hd Lmax 

are selected so that the difference between x^/DOF and 1 is the smallest. 

We get the results as follows. For defined by (@), (ijv = 0.28(1), A 
= 0.39(2), xVDOF = 1.21, and Q = 0.27. For m^' defined by (D, /S/u = 
0.313(9), A = 0.48(2), xVDOF = 1.22, and Q = 0.26. For defined by 
(§, j/u = 1.96(2), A = 0.040(3), xVDOF = 0.98, and Q = 0.47. For x^' 
defined by (D, 7/z/ = 2.37(2), A = 0.058(4), xVDOF = 1.21, and Q = 0.28. 
-^min = 36 and L^ax = 60 for m\ m)' , and x^' . Lmin = 30 and Lmax = 54 for 

The values of the critical exponents of w)' and x^' satisfy the scaling 
relation, 2/3 + 7 = dv where d is the dimensions, within errors. Thus we 
adopt them as the reliable estimates. Figure 1 shows the size dependence of 
w)' at z = Zc- The solid line indicates 0.48L~°'^^'^. Figure 2 shows the size 
dependence of x^' aX z = z^. The solid line indicates 0.058L^'^''. Figure 3 
shows the size dependence of U defined by (||) at z = z^- We cannot extract 
a reliable estimate in L = +00 since we do not know the correction to scaling 
of this system. It will have a value between 0.54 and 2/3. 
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4 Summary 



We carry out the Monte Carlo simulations of the hard-sphere lattice gas on 
the simple cubic lattice with nearest neighbour exclusion under fully periodic 
boundary conditions. Using the finite-size scaling, we obtain the critical 
exponents, P/v = 0.313(9) and 7/z/ = 2.37(2). They satisfy the scaling 
relation, 2/? + 7 = du. 

The critical exponents of the three-dimensional Ising model have been 
estimated, /?/z/ = 0.518(7) and 7/z/ = 1.9828(57) fS^]. These values do not 
agree with our results of the hard-sphere lattice gas. The value of the fourth- 
order cumulant of the Ising model on the simple cubic lattice under fully 
periodic boudary conditions is as follows. 0.44(2) [^], 0.4677(23) |3^, and 
0.4637(13) We see that these values are smaller than those of the hard- 
sphere lattice gas as seen in figure 3. While there is evidence that the hard- 
sphere lattice gas on the square lattice belongs to the Ising universality class 
in two dimensions as described in section |1|, it does not seem that the hard- 
sphere lattice gas on the simple cubic lattice falls into the universality class 
of the three-dimensional Ising model. 
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Figure captions 

Figure 1 Size dependence of tlie staggered density, m^', defined by (|^) at 
z = Zc- The solid line indicates 0.48L^°'^^^. Errors are less than the 
symbol size. 

Figure 2 Size dependence of the staggered compressibility, x^', defined by (|^) 
at z = Zc- The solid line indicates O.OSSL^'^''. Errors are less than the 
symbol size. 

Figure 3 Size dependence of the fourth-order cumulant, U, defined by (|^) 
at z = Zq . 
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Figure 3 
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